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^ Abstract 
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Evolution of the universe with modified holographic Ricci dark energy model is consid- 
Ci ered. Dependency of the equation of state parameter and deceleration parameter on the 

redshift and model parameters are obtained. It is shown that the density evolution of both 
^-H the non-relativistic matter and dark energy are same until recent times. The evolutionary 

^ trajectories of the model for different model parameters are obtained in the statefinder 

planes, r — s and r — q planes. The present statefinder parameters are obtained for different 
00 model parameter values, using that the model is differentiated from other standard models 

^ like ACDM model etc. We have also shown that the evolutionary trajectories are depend- 

ing on the model parameters, and at past times the dark energy is behaving like cold dark 
^ matter, with equation of state equal to zero. 

^ Keywords: Dark energy. Holographic model, Statefinder diagnostic, Cosmological evolu- 

tion. 

c3 PACS numbers: 98.80.Cq, 98.65.Dx 



1 Introduction 



Observations of distant type la supernovae (SNIa) and cosmic microwave background anisotropy 
have shown that the present universe is accelerating [1]. This expansion may be driven by a 
component with negative pressure, called dark energy. The simplest model of dark energy is the 
cosmological constant A which can fit the observations in a fair way [21 13], whose equation of state 
is = ~1- during the evolution of the universe. However there are two serious problems with 
cosmological constant model, namely the fine tuning and the cosmic coincidence To solve 
these problems different dynamic dark energy models have been proposed, with varying equation 
of state during the expansion of the universe. Holographic dark energy (HDE) is one among 
them [5|, El [7] . HDE is constructed based on the holographic principle, that in quantum gravity, 
the entropy of a system scales not with its volume but with its surface area L^, analogically 
the cosmological constant in Einstein's theory also is inverse of some length squared. It was 
shown that [5] in effective quantum field theory, the zero point energy of the system with size L 
should not exceed the mass of a black hole with the same size, thus Lr'pi^ < LMp, where pA is 
the quantum zero-point energy and Mp = 1/ a/SttG, is the reduced Plank mass. This inequality 
relation implies a link between the ultraviolet (UV) cut-off, defined through pA and the infrared 
(IR) cut-off encoded in the scale L. In the context of cosmology one can take the dark energy 
density of the universe px as the same as the vacuum energy, i.e. = Pa- The largest IR cut-off 
L is chosen by saturating the inequality, so that the holographic energy density can be written 
as 

p, = Sc^M^L-^ (1) 

where c is numerical constant. In the current literature, the IR cut-off has been taken as the 
Hubble horizon [HI C], particle horizon and event horizon [^ or some generalized IR cut off 
[El [HI HO] . The HDE models with Hubble horizon or particle horizon as the IR cut-off, cannot 
lead to the current accelerated expansion |6] of the universe. When the event horizon is taken 
as the length scale, the model is suffered from the following disadvantage. Future event horizon 
is a global concept of space-time. On the other hand density of dark energy is a local quantity. 
So the relation between them will pose challenges to the concept of causality. These leads to the 
introduction new HDE, where the length scale is given by the average radius of the Ricci scalar 
curvature, 

The holographic Ricci dark energy model introduced by Granda and Oliveros [11] based 
on the space-time scalar curvature, is fairly good in fitting with the observational data. This 
model have the following advantages. First, the fine tuning problem can be avoided in this 
model. Moreover, the presence of event horizon is not presumed in this model, so that the 
causality problem can be avoided. The coincidence problem can also be solved effectively in this 
model. Recently a modified form of Ricci dark energy was studied [12] in connection with the 
dark matter interaction, and analyses the model using Om diagnostic. In this paper we have 
considered the evolution of the universe in Modified Holographic Ricci Dark Energy (MHRDE) 
model and obtain the statefinder parameters to discriminate this model with other standard dark 
energy models. 

Statefinder parameters is a sensitive and diagnostic tool used to discriminate various dark 
energy models. The Hubble parameter H and deceleration parameter q alone cannot discriminate 
various dark energy models because of the degeneracy on these parameters. Hence Sahni et al. 
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[13] introduces a set of parameters {r, s} called statefinder parameters, defined as, 



''~aH^' '~ 3{q-ntotai)/2' 

where a is the scale factor of the expanding universe and fltotai is the total energy density 
containing dark energy, energy corresponds to curvature and also matter (we are neglecting the 
radiation part in our analysis). In general statefinder parameter is a geometrical diagnostic such 
that it depends upon the expansion factor and hence on the metric describing space-time. The 
r — s plot of dark energy models can help to differentiate and discriminate various models. For 
the well known ACDM model, the r — s trajectory is corresponds to fixed point, with r = 1 
and s = [13j. The cosmological behavior of various dark models including holographic dark 
energy model, were studied and differentiated in the recent literature using statefinder parameters 

[HIISIIISIE!. 

The paper is organized as follows. In section 2, we have studied the cosmological behavior 
of the MHRDE model and in section 3 we have considered the statefinder diagnostic analysis 
followed by the conclusions in section 4. 



2 The MHRDE model 



The universe is described by the Friedmann-Robertson- Walker metric given by 

ds^ = -dt" + a{tf ( + r^de^ + sin^ OdA , (3) 
\1 — kr^ J 

where (r, 6, 0) are the co-moving coordinates, k is the curvature parameter with values, k = 
1, 0, —1 for closed, fiat and open universes respectively and a{t) is the scale factor, with oq = 1, 
is taken as its present value. The Friedmann equation describing the evolution of the universe is 

i 

where we have taken 8ttG = 1, the summation includes the energy densities of non-relativistic 
matter and dark energy, i.e. J2i Pi = Pm + Px ■ The modified holographic Ricci dark energy can 
be expressed by taking the IR cutoff with the modified Ricci radius in terms of H and as 

HDIE] 

P.^^(h^'4h''), (5) 



a- (3 \ 2 

where H is the time derivative of the Hubble parameter, a and /3 are free constants, the model 
parameters. Chimento et. al. studied this type of dark energy in interaction with the dark 
matter with Chaplygin gas [I2], but our analysis is mainly concentrated on the cosmological 
evolution of MHRDE and analyses with statefinder diagnostic. Substituting the dark energy 
density as the MHRDE in the Friedmann equation, and changing the variable form cosmic time 
t to a; = In a we get 

ff=+*=^^ + _L^f! + ^i/= (6) 

3 3[a — p) dx a — p 
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Introducing the normalized Hubble parameter as h = H/ Ho and Qk 
Hubble parameter for x = 0, the above equation become, 



-2x 



'3x 



+ 



+ 



a 



3(a — (3) dx a — (3 



-k/HQ , where Hq is the 



(7) 



where Qmo = Pmo/3-f^o current density parameter of non-relativistic matter ( we will take 

0.27 as its values for our analysis throughout.) with current density p^o and flko is the present 
relative density parameter of the curvature. We will consider only flat universe, where Qko = in 
our further analysis. Solving the first order differential equation ([T]) we obtain the dimensionless 
Hubble parameter h as. 



''mo 



■3x 



1-/3 
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(a - /3)a, 



mO 



/3-1 



+ 1 



-3/3x 



Comparing this with the standard Friedman equation, the dark energy density can be identified 
as 



a 



1 



1-/3 



-3x 



+ 



{a - I3)n 



mO 



+ 1 



-3l3x 



(9) 



This shows that similar to the result obtained in references [H [15] for Ricci dark energy, the 
MHRDE density has one part which evolves like non-relativistic matter (~ e~^^) and the other 
part is slowly increasing with the decrease in redshift. The pressure corresponding the dark 
energy can be calculated as. 



Px 



-fix 



1 dVtx 
3 dx 



[{a-(3)n^o + (3-l\ e-'^" 



(10) 



Form the conservation equation, we can obtain the corresponding equation of state parameter 
for the fiat universe, using equation ^ as. 



-1 



IdhiVt^ 
3 dx 



■1 + 



{a - l)^]^o + /3 [(1 - /3) - (g - [3)n^o] e^d-/^)- 
(a - l)^]^o + [(1 -13) -{a- /3)^]„o] e^d-/^)- 



(11) 



This equation of state imphes the possibility of transit form cj^; > — 1 to a;^; < — 1, corresponds to 
the phantom model [TU] for suitable model parameter values. Recent observational evidences 
shows that the dark energy equation of state parameter can crosses the value -1 [20j. In a universe 
dominated with MHRDE, where the contribution from the non-relativistic matter behavior term 
is negligible in the dark energy density, the equation state parameter become. 



-l + (3 
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So if /3 is less than zero, the equation of state can crosses the phantom divide. In the far future 
of the universe, when redshift ^ — >■ — 1 also, the equation of state parameter reduces to the form 



given in equation (12). So the behavior of the dark energy is depending strongly on the model 
parameter [3. 



We have plotted the evolution of the equation of state parameter of MHRDE with redshift 
in figure [l| using the best fit values of the model parameters a and (3 as [12J (a, (3) = (1.01, 0.15) 
and r2mo=0.27. The evolution of Ux of dark energy shows that in the remote past of the universe, 
that is at large redshift, the equation of state parameter is near zero, implies that the dark energy 
behaves like the cold dark matter in the remote past. The plot also shows that at far future of 
the universe as z ^ —1, the equation of state parameter approaches a saturation value. The 
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Figure 1: Evolution of equation of state parameter oOx with redshift z for the best fit values 
a = 1.01 and (3 = 0.15. 
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Figure 2: Evolution of the equation of state parameter for other values of a and (3 

present value of the equation of state parameter according to this plot is negative, and is around 
ujx = -0.7. 



For other values of the model parameters [12j as (a,/3) = (1-01, 0.05), (1.01,0.1), 

(1.2,-0.05), (1.2,-0.1), (4/3,0.1) and (4/3, -0.1) the behavior of the equation of state parameter is 
given figure [2j For a given value of a the saturation value of Ux in the future universe decreases 
as 1/3 1 increases. 

Figure [2] shows that irrespective of the values of the parameters {a, (3), the equation of state 
parameter is negative at present times implies that the present universe is accelerating, and also 
in the remote past at high redshift Ux — 0, indicate that MHRDE behaves like cold dark matter 
in the past stages of the universe. For negative values of (3 the equation of state parameter 
crosses -1, in that case it can be classified as quintom [22] dark energy and for the case Ux < —I 
the universe will evolve into a phantom energy dominated epoch [21] . 

In figure |3} we have shown a comparison of the evolution of non-relativistic matter density 
and MHRDE density in logarithmic scale. Here we have neglected phase transitions, transitions 
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Figure 3: Evolution of non-relativistic matter density and MHRDE density in log. scale 

from non-relativistic to relativistic particles at high temperatures and new degrees of freedom 
etc. It is expected that these would not make much qualitative difference in the result. The plot 
shows that in the present model the densities of non-relativistic matter and dark energy were 
comparable with each other in the past universe that is at high redshift. The acceleration began 
at low redshifts, which solves the coincidence problem. 



The deceleration parameter q for the MHRDE model can obtained using the relation 



This equation can be expressed in terms of the dimensionless Hubble parameter h as 

1 dh^ 



2h^ dx 



1 



(13) 



(14) 



Using equation ([s] ) the above equation can be written as 
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(15) 



This equation shows the dependence of the deceleration parameter on the model parameters 
a and /3. As an approximation, if we neglect the contribution form the first terms in both 
numerator and denominator (since they are negligibly small) the deceleration parameter will 
become q = (3/3 — 2)/2. Which shows, as /3 increases form from zero, the parameter q increases 
form -1, that is the universe enter the acceleration phase at successively later times. In figure |4] 
and figure [5] we have plotted the evolution of the deceleration parameter with redshift. Figure |4] 
is for the best fitting model parameters q;=1.01, /3=0.15 and figure |5] is for the remaining model 
parameter values. The plots shows that at large redshift, the deceleration parameter approaches 
0.5. The universe is entering the acceleration in the recent past at 2; < 1. The plot also shows that 
as a increase, the entry to the accelerating phase is occurring at relatively lower values of redshift, 
that is the universe entering the accelerating phase at relatively later times as the parameter a 
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Figure 4: Evolution of deceleration parameter q for the best fit model parameters a=1.01 and 
/3=0.15 
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Figure 5: Evolution of deceleration parameter for other values of the model parameter. 



increases. The transition of the universe from deceleration to the accelerating phase is occurred 
at the the redshift Zt = 0.76, for the best fit model parameters . For comparison the combined 
analysis of SNe+CMB data with ACDM model gives the range Zt{ACDM) = 0.50 - 0.73 [201 El]. 
For taking consideration of the entire model parametric range, the transition to the accelerating 
phase can be obtained, as in figure [s] as Zt{MHRDE)=0. 50 - 0.76. The comparison of the two 
ranges shows that in the MHRDE model the universe entering the accelerating expansion phase 
earlier than in the ACDM model. The present value of the deceleration parameter for the best 
fit model parameters q;=1.01, /3=0.15 is go = —0.45 as from figure |4} 



3 Statefinder diagnostic 



We have calculated the statefinder parameters r and s, as defined earlier in equation ([2]). 
Statefinder parameters can provide us with a diagnosis which should unambiguously probe the 
properties of various classes of dark energy models. Equation ^ for r and s can be rewrite in 
terms of h"^ as, 

1 d^h^ 3 dh^ 
r = \ hi (16) 
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Figure 6: Evolutionary trajectory in the r — s plane for MHRDE model for the best fit values of 
the model parameters a=1.01, (3=0.15. The black spot on the top right corner corresponds to 
r = 1, s = the ACDM model. The today's point corresponds to r=0.59, s=0.15 
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(17) 



On substituting the relation for h"^ from equation (jSj), the above equations for a flat universe (in 
which = 0) become 



9/3(/3 



(a-/3)Q^o 
/3-1 



+ 1 e 



-3l3x 
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and s is become, 



/3(/3-l) 
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(a-/3)nmO 



1 e 



(a-/3)SlmO 
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+ 1 



(1 - /3)e-3/3^ 
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(19) 



From equations (18) and (19), it is evident that r = l, s = 0if/3 = and no matter what value 



a is, and this point in the r — s plane is corresponds to the ACDM model. This point is a very 
fixed point, thus statefinder diagnostic fails to discriminate between ACDM model and MHRDE 
model for the model parameter value (3 = 0. Since s is a constant for flat universe in this model, 
the trajectory in the r — s plane is a vertical segment, with constant s during the evolution of the 
universe, while r is monotonically decreasing form 1, if /3 is positive and monotonically increasing 
if P is assuming negative values. For a simple understanding, let us assume that Qm contribution 



is negligible small, when dark energy is dominating, then the equation (18) reduces to 

,-1 I W-1) 



(20) 



In this case for /3=0.05, 0.1 and 0.15 the corresponding values of r are 0.79, 0.60, 0.43 respectively. 
Bur when (3 assumes the negative values -0.05 and -0.10, the corresponding values of r become 
1.24 and 1.5 respectively. So at the outset the MHRDE model gives a r — s trajectory, as r 



starting form 1 and due to evolution of the universe the r will decreases to 1 



9/3(1-/3) 



if (3 is 



positive and increases to 1 



9/3(/3-l) 
2 ' 



if (3 is negative. 



The r — s evolutionary trajectory in the MHRDE model in fiat universe for the best fit model 
parameters a=1.01 and (3=0.15, is given in figure [6j In this plot as the universe expands, the 
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trajectory in the r — s plane starts form left to right. The standard ACDM model is corresponds 
to r = 1, s = is denoted. In this model the parameter r first decrease very slowly with s, then 
after around s=0 r decreases steeply. The today's value of the statefinder parameter (ro=0.59, 
So=0.15) is denoted in the plot. 
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Figure 7: the first plot is for a,/3 = 1.01,0.10 and the second plot is for a,f3 = 1.01,0.05 The 
black spot on the top right corner corresponds to ACDM model, the present state of the evolution 
is denoted as today's point. 

For other model parameters, the r — s plots are given in figure [7j These plots also shows 
the same behavior of figure |6| but the separation between ACDM model and MHRDE model 
in the r — s plane is increasing as /3 increases. The respective todays universe corresponds 
ro. So = 0.71, 0.1 and tq, Sq = 0.85, 0.05. 



For negative values of /3, the evolutionary characteristics is plotted in figure |8] for model 
parameters a, /3 = 1, 2, —0.10; 4/3, —0.10. Here also the evolution in the r — s plane is from left 
to right. In this case the behavior is different form that for the positive /3 value, in the sense that 
as s increases The r is increasing to vales greater than one. The increase is very slowly at first 
then increases steeply as the universe evolves. The today's value in these cases are ro=1.325, 
so=-0.10 when (a,/3)= (1.2,-0.10) and ro=1.321, so=-0.10 for («,/3) = (4/3, -0.10) respectively. 
The difference between MHRDE model for these model parameters and ACDM can be noted. 

The statefinder diagnostic can discriminate this model with other models. As example, for 
the quintessence model the r — s trajectory is lying in the region s > 0, r < 1 and for Chaplygin 
gas the trajectory is in the region s < 0,r > 1. Holographic dark energy with the future event 
horizon as IR cutoff, starts its evolution form s = 2/3, r = 1 and ends on at ACDM model fixed 
point in the future [Ml l27j . 

In order to confirm the r — s behavior of MHRDE model, we have plotted the behavior in 
r — z plane, in figure |9j For MHRDE model, the r value is commencing from 1 irrespective of 
the values of a and f3 at remote past and as the universe evolves, r is decreasing, if /3 is positive 
and increasing, if (3 is negative. 



We have studied the evolutionary behavior in the r 



q plane also. For positive value of {3 
the plot is shown in figure 10 for the standard values a = 1.01, /3 = 0.15 . The figure shows 
that both ACDM model and MHRDE model commence evolving from the same point in the 
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Figure 8: r — s plots for model parameters a, /3 = 1.2,-0.10; 4/3,-0.10. the arrow in lower 
left corner of the panel shows the evolution towards ACDM model. The present position of the 
evolution is denoted as today's point. 
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Figure 9: r — z plot, all model parameters. Shows that for positive values of /3 r is decreasing 
form 1, but for negative values of (3 the value of r is increases form 1. 
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Figure 10: Evolutionary trajectory in the statefinder r — q plane with a — 1.01 and p — 0.15. 

The solid line represents the MHRDE model, and the dashed line the ACDM (denoted as LCDM 
model in the plot) as comparison. Location of today's point is (0.59, -0.45). 
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Figure 11: Evolutionary trajectory in the r — q plane with a=1.2, (3= -0.10. The present position 
is denoted. The dashed shows the evolution of ACDM (denoted as LCDM model in the plot) 
model from right to left. 

past corresponds to r = 1, g = 0.5, which corresponds to a matter dominated SCDM universe. 
In ACDM model the trajectory will end their evolution at g = — l,r = 1 which corresponds to 
de Sitter model, while in MHRDE model the behavior is different from this. The statefinder 
trajectory in holographic dark energy model with future event horizon has the same starting 
point and the same end point as ACDM model [251 ESj. Thus MHRDE model is also different 
form holographic dark energy with event horizon form the statefinder viewpoint. 



For negative values (3 the plot is as given in figure 11 The evolution of the trajectory 



is starting from left to right. Note that the r value is at the increase from one as the universe 
evolves. It is evident from the plot that the present position of the model corresponds to ro=1.325 
and go = -0.63. 



4 Conclusions 



We have studied the modified holographic Ricci dark energy (MHRDE) in fiat universe, where 
the IR cutoff is given by the modified Ricci scalar, and the dark energy become px = 2 (if + 
3aH^ /2)/{a — (3) where a and (3 are model parameters. We have calculated the relevant cos- 
mological parameters and their evolution and also analyzed the model form the statefinder view 
point for discriminating it from other models. The importance of the model is that it depends 
on the local quantities and thus avoids the causality problem. 

The density of MHRDE is comparable with the non-relativistic matter at high redshift as 
shown in figure |3] and began to dominate at low redshifts, thus the model is free from the 
coincidence problem. 

The evolution of equation of state parameter is studied. The equation of state parameter 
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is nearly zero at high redshift, imphes that in the past universe MHRDE behaves hke cold dark 
matter. Further evolution of equation of state is strongly depending on the model parameter /3. 
If the P parameter is positive the equation of state is greater than -1. For negative values of f3, 
the equation of state cross the phantom divide < —1. 

In this model the deceleration parameter starts form around 0.5 at the early times and 
and starts to become negative when the redshift z < 1. . In general we have found that in 
MHRDE model the universe entering the accelerating phase at times earlier (for allowed range 
of parameters a and /3), than in the ACDM model. But in particular as the model parameter a 
increases, the universe enter the accelerating phase at relatively later times. 

We have applied the statefinder diagnostic to the MHRDE and plot the trajectories in the 
r — s and r — q planes. The statefinder diagnostic is a crucial tool for discriminating different dark 
energy models. The statefinder trajectories are depending on the model parameters a and /3. For 
positive values of (3 the r values will decreases from one and for negative /3 the r will increases 
form one as the universe evolves. The values of a and /3 are constrained using observational data 
in reference [12], the best fit value is a=1.01, /3=0.15. The present value of (r, s) can be viewed 
as a discriminator for testing different dark energy models. For the ACDM model statefinder is 
a fixed point r=l, s=0. For positive values f3 parameter the r — s and r — q plots of MHRDE 
shows that, the evolutionary trajectories starts form r = 1 and q = 0.5, in the past universe 
(for the best fit model parameters), which reveals that the MHRDE is behaving like cold dark 
matter in the past. The further evolution of MHRDE in the r — s plane shows that the present 
position of MHRDE model in the r — s plane for the best fit parameter is ro=0.59, so=0.15 
and in the r — q plane is ro=0.59, go=-0.45. The difference between the MHRDE and ACDM 
models is in the evolution of the equation of state parameter, which is -1 in the ACDM model 
and a time-dependent variable in MHRDE model. A further comparison can be made with the 
new HDE model [H], which gives the present values rQ^HDE) = 1.357, Sq{HDE) = —0.102 and 
ro{HDE) = 1.357, qo{HDE) = -0.590. So in the r - s plane the distance of the MHRDE model 
form the ACDM fixed point is slightly larger compared to the new HDE model for positive values 
of P parameter. However in the case of MHRDE model the starting point in r — s plane and 
r — q plane is (r = 1, s = and r = l,q = 0.5) is same as that in the ACDM model. 

For negative values of the (3 the r — s trajectory we have plotted is different compared to 
that of positive /3 values. For negative /3 values the r value can attains values greater than one 
as s increases. The present status of the evolution in the r — s plane is tq =1.325, so=-0.10 for 
model parameters a=1.2, /3=-0.10 and ro=1.321, so=-0.10. The r — q for a=1.2 and /3=-0.10 
shows that the present state of the MHRDE model is corresponds to tq = 1.325 and qo= - 0.63. 
These values shows that the MHRDE model is different form ACDM model for the present time 
when /3 parameter is negative also. But compared the new HDE model, the present MHRDE 
model doesn't show much deviation, shows that for negative values the behavior of MHRDE 
model is almost similar to new HDE model. Irrespective of whether /3 is positive or negative 
the MHRDE model is commence to evolve from SCDM model. When /3 is positive 1 is the 
maximum value of r, on the other hand when /3 is negative 1 is the minimum value of r. However 
the exact discrimination of the dark energy models is possible only if we can obtain the present 
r — s values in a model independent way form the observational data. It is expected that the 
future high-precision SNAP-type observations can lead to the present statefinder parameters, 
which could help us to find the right dark energy models. 
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